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Abstract 

Due to the smallncss of the lepton Yukawa couplings, higher dimensional operators can 
give a significant contribution to the lepton masses. In this case, the lepton mass matrix 
and the matrix of lepton-Higgs couplings are misaligned leading to lepton flavor violation 
(LFV) mediated by the Standard Model Higgs boson. We derive model-independent bounds 
on the Higgs flavor violating couplings and quantify LFV in decays of leptons and electric 
dipole moments for a class of lepton-Higgs operators contributing to lepton masses. We find 
significant Higgs-mediated LFV effects at both 1-loop and 2-loop levels. 



1 Introduction 



The flavor puzzle of the Standard Model remains one of the outstanding problems in particle 
physics. Masses of elementary fermions range over many orders of magnitude, forming a pattern 
which cannot be explained within the Standard Model. Also, it is not known whether these 
masses are generated by renormalizable terms in the Lagrangian or higher order operators also 
contribute pjj [2]. This issue will be studied at the LHC by means of the Higgs couplings 
measurements [3], which will give us at least partial answers. Complementary information 
about possible effects of non-renormalizable operators is provided by observables sensitive to 
flavor violation. 

In this work, we focus on the lepton sector due to its particular sensitivity to flavor violation. 
We will consider a class of higher dimensional operators involving the Higgs field which affect 
the lepton masses. Although such operators are suppressed by a "new physics" scale M, they 
can still give a significant contribution to the lepton masses due to the smallness of the lepton 
Yukawa couplings. In particular, we will focus on operators of the type 

- AC = Y.f H L Li l Rj + Y£ )H ^ H L Li l Rj + ... + h.c. , (1) 

as well as higher order operators. Here and Y-j are a priori independent flavor matrices. 
Since the resulting lepton mass matrix and the matrix of the Higgs couplings are in general 
misaligned in flavor space, the presence of higher order operators leads to flavor violation me- 
diated by the Standard Model Higgs boson. The rotation to the mass eigenstate basis may 
involve large angles, as hinted by the neutrino sector, which amplifies the effect. For example, 
the processes involving light generations such as [i — > ej can be dominated by loops involving 
the t lepton. The existing bounds on LFV observables then place strict constraints on the scale 
of new physics and/or the type of admissible Yukawa textures. 

Some aspects of lepton flavor violation in Higgs interactions have been considered before (for 
early studies, see [U El E]), although a systematic study is lacking. For example, LFV Higgs 
decays h — > klj were analyzed in 0IH]- The Higgs-mediated decay r — >• fx'j was considered in 
[9], where it was found that this mode is not particularly constraining. Effects of higher order 
-operators on neutrino masses and LFV in a 2 Higgs doublet model were studied in |10j . 
Finally, Ref. [11] summarizes Higgs-induced quark FCNC effects in extended Higgs models. 

The paper is organized as follows. In Section 2, we derive model-independent bounds on the 
flavor violating couplings of the Higgs boson. In Section 3.1, we apply these bounds to the case 
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of dimension 6 operators and obtain constraints on the scale of new physics and on the rotation 
angles. In Section 3.2, we study the possibility that the lepton mass hierarchy is created entirely 
by non-renormalizable operators. In Section 4, we present our conclusions. 

2 Bounds on the Higgs couplings 

The relevant Lagrangian describing interactions of the physical Higgs boson h with leptons is 
given by 

AC = -^L h kPRlj + h.c. , (2) 
v2 

where Pr = (1 + 75)/2. This interaction induces the following flavor changing dipole operator 
(Fig.0, 

£ cffl = eL i:j ka^F„ v P L lj + h.c. , (3) 
where = |[ 7 M ,7l, Pl = (1 - 7s)A and 

Li i = aA 2 2 "ir In — . 4 

Here we have kept only the leading in m T /mh r-lepton contribution in the loop since the others 
are negligible for our purposes. The corresponding branching ratio is given by 

BR(lj l a ) = BR(lj -> kvjVi) x l ^ ° (|L y | 2 + (Ljil 2 ) . (5) 

The Higgs interactions also induce flavor-diagonal dipole operators at one loop: 

£ eff2 = e ReLjj ka^F^Ji - ie ImLu ka^F^^k + h.c. . (6) 

These are constrained by the charged lepton anomalous magnetic moments and electric dipole 
moments, 

|<5aj| = 4mj|ReLjj| , 
\di\ = 2e|lmLjj| . (7) 



Finally, there are tree-level processes lj — > Ulkl^ induced by ([2j). Their branching ratio is 
given by 

BR(lj -> kl k l+) = BR(lj -> hujPi) x 2 ( 5 4 6 ~ 2 ^ 4 lto| 2 (|^f + |%i| 2 ) • (8) 

_P h 

The resulting bounds on the Higgs couplings are presented in Table 1. Most experimental 
constraints are taken from Particle Data Group [12], while the recent bound on BR(/i — > e-y) is 
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Figure 1: Leading Higgs-mediated contributions to the lepton dipole operators and 3-body 
decays. 




Figure 2: Leading 2-loop Higgs-mediated contribution to the lepton dipole operators. 



from [13J and the interpretation of 5a e is due to [H]. The third column shows the combination of 
Higgs couplings constrained by a particular observable, while the fourth one shows representative 
bounds on the couplings under the assumptions j/y = yji and yu = rm/v, as in the Standard 
ModelJE 

2.1 2— loop contributions 

Two loop diagrams of certain types can be important as they are suppressed by only one power 
of small Yukawa couplings [4]. In what follows, we include the leading Barr-Zee type diagrams 
[15] with the top quark in the loop (Fig. [2]) [16] . Since the coupling of the top quark to the 
Higgs is essentially unaffected by the higher dimensional operators, we will assume the SM value 
yu = 1. Then [16], 

BR(Zj -> la) = BR(lj -> UvjVi) x f( z ) (\Vij\ 2 + 1%'il 2 ) » ( 9 ) 



where d = 174 GeV and 



„ . . 1 Z" 1 , 1 - 2»(1 - a;) x(l-x) 

f(z) = -z / dx \ >- In ^ >- , 10 

2 ./ x(l — x) — z z 



1 Here we have not included the bound from fi to e conversion on nuclei since it depends on the Higgs couplings 
to light quarks, which are uncertain in our framework. In any case, for the SM Higgs-quark couplings, the bound 
is of order \/\yi2\ 2 + |j/2i| 2 < 5 x 1CP 4 for m h — 200 GeV, which is superceded by the 2-loop Barr-Zee bound to 
be studied below. 



observable 



present limit 



constraint 



y3i2/23| 2 + |y32yi3| 2 J < 7 x 

\ 1/4 

2/33I 2 (Irf + li/d 2 )) < 5 x io- 2 
y 33 | 2 (|y 3 i| 2 + |yi3l 2 )) 1/4 <5xio- 2 
yn\ 2 (|y 2 i| 2 + |yi 2 | 2 )) 1/4 <2xio- 3 
|y 22 | 2 (|y 23 | 2 + |y 32 | 2 )j < 4 x io~ 2 

\ 1/4 

\m? (Irf + bil 2 )) < 4 x io- 2 
2/22I 2 (|yi3l 2 + |y 3 i| 2 )) 1/4 <4x io- 2 



constraint for 

Vij = Vji , yu = m/v 



y/\yi3m\ < 6 x io~ 4 

\y 23 \ < 2 x HT 1 
|yi 3 | < 2 x HT 1 
1 2/12 1 < 1 
1 2/23 1 < 1-7 
1 2/13 1 < O(10 2 ) 
1 2/13 1 < 1.7 



BR (n -»• e 7 ) 
BR (r -»• /i7) 



BR (r 
BR (ju 
BR (r 
BR (r 
BR (r 



e7) 
eee) 

eee) 
enn) 



2.4 x 10~ 12 
4.4 x 10~ 8 
3.3 x IO" 8 

1.0 x IO" 12 

2.1 x 10~ 8 
2.7 x 10~ 8 
2.7 x 10~ 8 



c/ e (e-cm) 
(e-cm) 

(5a e 
<5a„ 



1.1 x 10~ 27 
3.7 x 10~ 19 
2.3 x IO -11 
40 x IO" 10 



|lm(y 3 i2/i3)| < 2 x IO" 4 
|lm(j/ 3 22/23)| < 4.1 



Im(y 2 3 ) < 2 x IO" 4 



Im(y 2 



23 , 



< 4.1 



v ,Re(y 3 iyi3)| < 0.14 
^|Re(j/ 32 2/ 23 )| <0-13 



Re(y 2 3 )| <0.14 



Re(y 2 3 )| <0.13 



Table 1: Current experimental limits on flavor and CP violating observables in the lepton 
sector, and the corresponding constraints on the Higgs couplings. The displayed bounds on y^ 
correspond to rrih = 200 GeV; for other Higgs masses they are to be multiplied by m/j/(200 
GeV). 



observable 


constraint 


BRO -»• erf) 
BR (r -»• nj) 

BR (r -»• e7) 

d e 

dp, 
5a e 
5 a p 


(l2/i 2 | 2 + l2/ 2 i| 2 ) 1/2 <6xl0- 6 

/ \ I/ 2 

(jrf + \V32\ 2 ) <4xl(T 2 

/ \ I/ 2 

(l2/i3| 2 + |y 3 i| 2 ) <3xio- 2 
|lm yu\ < 6 x 10~ 7 
|lm y 22 | < O(10 2 ) 
Re y u < 0.27 
Rey 22 <0.21 



Table 2: Two loop constraints for = 200 GeV. (The bounds get tighter by a factor of 1.5 
for m h = 100 GeV.) 
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with z = mf/m?. The resulting constraints are presented in Table 2. Since f(z) varies from 
0.76 for rrih = 200 GeV to 1.14 for = 100 GeV, these bounds are not very sensitive to the 
Higgs mass. 

The flavor diagonal transitions can also be extracted from the amplitude computed in |16j . 
We find 

Ul = 24^3 /(*) Vti ■ (11) 
The EDMs and g — 2 then constrain the diagonal couplings as summarized in Table 2. We note 
that the bounds from BR(/x — > ej) and d e are particularly important^] 

3 Lepton flavor violation from Yukawa— type interactions 

3.1 Dimension— 6 operators 

The lepton sector is particularly sensitive to BSM flavor structures due to the smallness of 
the lepton masses and strong constraints on lepton flavor violation. Since the lepton Yukawa 
couplings in the SM can be as small as 10~ 5 , higher dimensional operators involving the Higgs 
field can give a significant contribution to the lepton masses. In this section we consider an 
effect of dimension 6 operators of this sort, 

-AC = H L Li l Rj (Vf + Y<V ^fj + h.c. , (12) 

which amounts to replacing the constant SM Yukawa couplings with Higgs-dependent ones, 

Here , are in general independent flavor matrices and M is the new physics scale. An 
immediate consequence of the above Lagrangian is that the SM Higgs boson mediates tree level 
flavor changing neutral currents. Indeed, the lepton mass matrix is given by 

Mv=v (y<P + Y<^, (14) 

whereas the matrix of couplings of the physical Higgs boson is 

where we have used the convention H° = v + h/y/2. Clearly, these matrices are in general 
misaligned in flavor space and the Higgs couplings in the mass eigenstate basis (cf. Eq. [2J can 
2 In 2 Higgs doublet models, Barr-Zee contributions often dominate at low tan/3 |17j . 
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change flavor. The latter are given by 

y = ulyU R , (16) 

where the unitary matrices Ul, Ur diagonalize the lepton mass matrix, M C//j=diag(m e , m^, m T ) 

To understand the constraints on the dim-6 operators, we vary the proportion between 
and Y^ ] v 2 /M 2 and scan over different Ul, Ur. We allow for large angle rotations, as hinted by 
the neutrino sector, while being agnostic about the neutrino mass matrix which presumably is 
provided by the seesaw mechanism. The resulting yij are then subject to the bounds described 
in the previous section. 

For the scan, we use the following procedure. The Yukawa textures are generated through 

Y = U L ~ di&g(m e ,mu,m T ) ut , (17) 
v 

by scanning over Ul, Ur. Modulo phase redefinitions of the lepton fields, Ul,r can be chosen as 

U L = V L , Ur = V r @, (18) 

with unitary Vl,r parametrized by three mixing angles and a single phase as in Refs. |12} [T8] , 
and being a diagonal phase matrix, 

e = diag(e^ 1 ,e^ 2 ,e^ 3 ) , (19) 

subject to the constraint <pi + <p2 + 4*3 = 0. Note that there are 4 reparametrization-invariant 
phases in the lepton mass matrix since 5 out of 9 original phases can be eliminated by the phase 
transformations of the left-handed and right-handed leptons. 



Having generated Yij (setting H — > v), we split it into two pieces according to Eq. 13 Clearly, 



if Ytf dominates, Mjj and 3^ are almost aligned, and the lepton FCNC are suppressed. The 
latter therefore set a bound on the allowed proportion of Y i jV 2 /M 2 in the lepton mass matrix. 
We scan over the following parameters: 

• 6 angles of Vl, Vr 

• 2 phases of Vl, Vr and 2 phases of G 

• 9 complex parameters 
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In what 



Note that, given Y^ and Yj Uj , the remaining piece Y^ v 2 /M 2 is determined by Eq. 
follows we present our results for two cases: Y^ is varied first in the range [0.5, 1] x Yy and 
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Figure 3: BR(/j — > e7) vs |d e | (left) and BR(/x — > e7) vs BR(/i — >• eee) (right) for arbitrary 
rotation angles and rrih = 200 GeV. Dimension-6 operators contribute up to 50% (top) and 10% 
(bottom) to Yij. The experimental limits are given by the straight lines. 

then in the range [0.9, 1] x This restricts the relative contribution of the dim-6 term to no 
more than 50% and 10%, respectively. 

Our results are presented in Fig. [3j We see that if the dim-6 operators are allowed to 
contribute as much as 50% to Yij, BR(/i — > e'y) and d e are overproduced by up to 2 orders of 
magnitude. On the other hand, if this contribution is below 10%, most points are allowed for 
rrih = 200 GeV. Thus, allowing for arbitrary rotation angles, we find an empirical constraint 

.„ , 200 GeV . . 

< 0.1 x . 20 



'« M 2 

This can be reinterpreted in terms of the bounds on the new physics scale M. For the two 



3 This range applies separately to the real and imaginary parts of ■ 
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limiting cases of similar Y^y and Ya, and order one Yjy , we get 

Y^~Y l3 =► M> 500 GeVx^^X, 

M>200 TeVx 2 ^X. (21) 

In the latter case, we used the most restrictive Yukawa couplings 0(1O -5 ) involving the electron. 
Let us note that we find observables other than BR(/i — > ej) and d e far less constraining and 
also confirm numerically that the diagrams with muons/electrons in the loop are unimportant. 
The latter statement is easy to understand. For the r contribution, the relevant Higgs vertices 
involve the r mass times a mixing parameter, e.g. 

2/13 ~ ei3 — , (22) 
v 

and, in addition, the mass insertion in the loop is m T . This enhances the r loop by orders of 
magnitude. The couplings involving only the first two generations are typically bounded by 
m n/ v ~ 10 -3 , unless there is a very strong mixing with the r. To account for the muon loop 
contribution to BR(/x — > ej) and d e , one can simply rescale the bounds of Table 1 by \Jm T /m^ 
and replace index 3 by 2. Keeping in mind that there are mixing angles appearing at the vertices 
so that the actual couplings are smaller than 10 -3 , one finds that these constraints are satisfied 
automatically. 

As the next step, we study constraints on the rotation angles parametrizing Ul,r, while 
allowing for arbitrary values of the four phases therein as well as a wide range of proportions 



between Y^ ] and Y^v 2 /M 2 in Eq. 13 Specifically, we choose Y^ in the range ±0.9 x Y^. 
The results are presented in Fig. |4j where we vary all the angles (in Ul and Ur) in the same 
range. If one allows for angles as large as 0.1, both BR(/j — > cy) and d e are overproduced, while 
the angles of order 0.03 are typically consistent with the constraints for = 200 GeV. More 
precisely, BR( / u — > ej) and d e are most sensitive to the 1-3 and 2-3 mixing angles #13, #23 in 
terms of the standard parametrization of unitary matrices [TH]. We then find 

0.3,0,3 < 3* 10- >< 5^ , (23) 

while 612 is allowed to be as large as 0(1). We note that d e is somewhat more restrictive than 
BR(/i ->■ ej) and requires 0i 3 , 623 < 10~ 2 for m h = 200 GeV. 



(0) 
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Figure 4: BR(^ — > ej) vs \d e \ for small rotation angles 9. Left: 9 < 0.1, right: 9 < 0.03. Here 
rrih = 200 GeV. The experimental limits are given by the straight lines. 

3.1.1 2 loop constraints 

The 2-loop contributions are very sensitive to the couplings in the (12)-block. In particular, 
if one assumed a "natural" value of yi2 to be of order ^/m e m^/v, BR(^ — > ej) would be 
overproduced by orders of magnitude, as seen from Table 2. The resulting constraints are very 
strong, as our numerical study demonstrates in Fig. [5j An analog of Eq. 20 in the (12)-sector 
reads 



a m 2 



< W-^Yijl (24) 



for nih = 200 GeV. The Higgs mass dependence is weaker than linear and the bound gets tighter 
by a factor 1.5 for mh = 100 GeV. Interpreting this constraint as is done in Eq. [2T| we get 

~ Yij M > 5 TeV , 
Yfp ~ 1 M > 2000 TeV . (25) 

Allowing for arbitrary proportions between dim-4 and dim-6 operators, we find the following 
constraints on the rotation angles 

012 < 10~ 3 , 0i3,23 < 1(T 2 . (26) 

Note that BR(/i — > cy) is also sensitive to #13,231 although not as much as to 9\2- Unlike in the 
1-loop case, a diagonal lepton Yukawa matrix does not solve all the problems: d e is sensitive to 
the phase of the flavor-diagonal coupling yn which does not vanish in the limit of zero rotation 
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Figure 5: BR(^ — > ej) vs \d e \ at 2 loops. Left: the dim-6 contribution is limited to 0.1%, while 
the rotation angles are arbitrary. Right: the rotation angles are bounded by 6±2 < 10~ 3 , #13,23 < 
10~ 2 and the CP phases 4> < 0.1, while the balance between dim-4 and dim-6 contributions is 
arbitrary. Here rrih = 200 GeV. The experimental limits are given by the straight lines. 

angles. We thus have a separate constraint on the CP phase of y\\: 

4> < W- 1 . (27) 
3.2 Higher dimensional operators and the mass hierarchy 

Let us now explore the possibility that the lepton mass hierarchy is entirely due to higher 
dimensional operators [IJ [2] . The Yukawa couplings are expanded as 



00 



n,,=0 V 7 

with order one lJ and M being a new physics scale. In most interesting cases which address 
the flavor problem, the coefficients 1 vanish up to a certain order mj. This can happen 
due to some symmetry (e.g. Froggatt-Nielsen type |19j ) of the UV completion of our effective 
theory, which may not be apparent at low energies. In this case, the mass hierarchy is generated 
by 

v 2 



and the Yukawa textures take the form 



" = W <<1 (29) 



Yij = c i3 , (30) 
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with order one Qj . Analogous textures in the quark sector were considered in [2] . It is interesting 
that the invariant measure of CP violation increases by many orders of magnitude compared to 
that in the Standard Model, which can be relevant to baryogenesis |20j . 

Various lepton textures of the above type can be generated as follows. The Yukawa matrix 



is represented in terms of the eigenvalues and the rotation matrices Ul,r as in Eq. 17 Scanning 
over Ul,r then produces viable lepton textures. Choosing for definitness e = 1/60 as motivated 
by the top-bottom quark mass hierarchy [2], one determines the exponents via 

riij = round(log e \Yij\) , (31) 

which also fixes the 0(1) coefficients c^-. Moving to the mass eigenstate basis, we obtain the 
lepton couplings of the physical Higgs ^ , 

Vii = (U L y ki (2n« + 1)Y« (U R )ij . (32) 

Note that, in this basis, yij are defined up to a phase (see [20] for details). Indeed, one can 
multiply the left- and right-handed leptons by diagonal phase matrices such that the mass 
terms remain the same. However, observables are invariant under this reparametrization since 
they involve combinations like y^iUa ° r absolute values of the couplings. 

For concreteness, let us restrict ourselves to Yukawa textures of the following type 



e 3 






e l 


e 2 


e n 




e n 


e 1 



7~ f e 2 e n . (33) 
\e m e n e 1 / 

Here we have chosen the diagonal entries such that they reproduce (roughly) the m e : : m T 
hierarchy. The integers I, m, n are subject to the LFV and CP constraints. For large /, m, n the 
Yukawa matrix is approximately diagonal and the constraints are satisfied. Numerically, we find 
the following approximate bound: 

I > 4 , m>2 , n>2. (34) 

For most choices of order one coefficients in the texture, the constraints are satisfied in this case. 
Note that the bound on I is the strongest one due to the Barr-Zee contributions to — > erf. 
This entails, in particular, that the lepton analog of the down quark texture studied in Ref. [2\ 
is strongly disfavored. 

Our numerical results for representative textures are shown in Fig. [6} We take I = 4, 
m = n = 2 for panel (a) and I = m = 2, n = 1 for panel (b). In this figure, the light 
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Figure 6: BR(/j — > ej) vs \d e \ for representative Yukawa textures, with m/j = 200 GeV. Light 
(orange) and dark (blue) dots correspond to one loop and one + two loop contributions, respec- 
tively. The experimental limits are given by the straight lines. 

(orange) and dark (blue) dots correspond to one loop and one loop plus two loop contributions, 
respectively. For case (a) most points satisfy the constraints, whereas in case (b) the bounds are 
grossly violated. We find that BR(/u — > e~f) is dominated by the Barr-Zee contributions, whereas 
d e receives comparable contributions at 1 and 2 loops. As expected, significant departures from 
the diagonal Yukawa form, especially in the (12)-block, are not allowed. 

Finally, let us note that although the constraints are strong, the remarkable sensitivity of 
/i — > &y and d e also implies good prospects for detection of Higgs-induced lepton flavor and CP 
violation in current and future experiments. 

4 Conclusion 

In this work, we have derived model-independent bounds on the lepton flavor violating cou- 
plings of the SM Higgs boson. Such flavor violation appears when higher dimensional operators 
contribute to the lepton masses. Scanning over various textures shows that these contributions 
should be limited to about 0.1% in the (12)-sector and no more than 10% in the (13) and 
(23)-sectors, otherwise BR(^ — > e^f) and d e are overproduced. Alternatively, if one allows for 
large contributions of the higher order operators, the Yukawa matrix must be diagonalizable by 
a small-angle rotation with 812 ~ O(10~ 3 ) and #13, #23 ~ 0(1O -2 ). In addition, the electron 
EDM sets a constraint on the relevant CP phase 4> < O(10 _1 ). 
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Further, we have studied the possibility that the lepton mass hierarchy is created entirely by 
non-renormalizable operators. Also in this case the LFV effects are significant. The preferred 
textures have small intergenerational mixing, especially in the (12)-sector, as required by the 
Barr-Zee contributions to fj, — >• ej. As a result, the lepton analog of the down quark texture 
studied in [2j is strongly disfavored. 

On the other hand, the remarkable sensitivity of LFV and CP violating observables to Higgs- 
induced effects implies good experimental prospects for detection of \x — > e^f and d e . 
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